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SUMMARY 

A procedure is presented for the design of two-dimensional 
channels for compressible nonviscous flow. The method requires 
boundary conditions consisting of the shape of one c han nel sur- 
face and the velocity distribution on that surface. The process 
consists of the step-by-step computation of an arbitrary number of 
streamlines within the channel. Two variations of the method are 
presented; One variation is based on an assumed constant stream- 
line curvature along an equipotential line, and the other is based 
on an assumed vortex -type variation. Tables for the use of both 
methods are presented for a range of Mich numbers and values of 
the ratio of specific heatB. 


INTRODUCTION 

Hi order to obtain satisfactory performance from gas turbines, 
blade-profile design yielding prescribed' velocity distributions is 
desirable. Proposed procedures for blade-profile design from 
specified surface velocities known as inverse procedures, such as 
those given in references 1 and 2, utilize conformal transformation 
or interference techniques. Such methods are limited in application 
by the assumption of incompressible flow and by the lengthy com- 
putations required. 

For closely spaced blades of hi$i solidity, such as those 
encountered in axial -flow turbines, the application of channel 
methods has been found useful. Stodola (reference 3) presents 
Flugel’s stream-filament technique for determining the velocities 
in a known channel by graphical integration. A simplified method 
for such analysis is presented in reference 4. 

A procedure based an stream-filament techniques for the 
design of two-dimensional channels for compressible nonviscous 
flow has been developed at the NACA Lewis laboratory and is 
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presented herein. The method requires boundary conditions con- 
sisting of the share of one channel surface and the velocity 
distribution on that surface. Although the method is not a 
completely inverse solution, it is useful in the design of blades 
for hi solidity cascades. The proposed procedure deals only 
■with the channel between adjacent blades and is inapplicable to 
the regions in the vicinity of the nose and the tail. A similar 
technique for incompressible flow is given by von Mises (refer- 
ence 5). The specification of the shape and the velocity dis- 
tribution on one channel surface constitutes an over-determination 
of desired conditions in that the specification of only a short 
continuous length of streamline and the velocities along it are 
needed to determine the flow in the complete field. The method of 
solution, however, is based on satisfying the conditions of 
continuity and irrotational flow at a series of points along the 
specified streamline and the fairing of a curve to represent the 
solution. The solution is thus an approximate one and the 
velocities in the resulting channel should be computed to determine 
the degree of satisfaction of the specified conditions. 

Further variation from specified conditions will result from 
the fact that the inlet flow required to satisfy theoretically the 
specified conditions in the channel is unknown. This variation will 
undoubtedly be small as experiments of high-solidity cascades have 
shown that the flow in the channel sections of cascades is relatively 
insensitive to the inlet conditions. 

In the application to cascade design of the method presented 
herein, specification of the boundary conditions as those for the 
suction surface of the blade is desirable. An arbitrary number of 
stream filaments is successively constructed starting from tne 
assigned surface, each carrying a specified fraction of the total 
weight flow. The streamline bounding the last stream filament is 
taken as the channel boundary and as the pressure surface of an 
adjacent blade*. The orthogonal spacing between a given streamline 
and the next streamline constructed from it is approximated by a 
spacing taken perpendicular to the given streamline. This spacing • 
is computed by integration of the flow equations with an assumed 
variation of streamline curvature alo n g an equipotential line. 

Two variations are considered herein $ (l) a constant curvature, 
and (2) a vortex-type variation. 

Weinig (reference 6) presents a solution to a similar problem 
by the graphical construction of an orthogonal flow net in which 
the spacings between streamlines, for equal potential increments 
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along -til© assigned streamline, are taken as equal to the potential 
spacings with corrections for compressibility. The streamlines 
are constructed for several potential increments an d an extra- 
polation yields the final flow net. 


THEORY CEF METHOD 

The problem may be stated as follows : Given the shape of one 

of the channel boundaries, the desired velocity distribution therecn, 
and the weight flow to be carried by the channel, find the. opposite 
boundary. In the case of a cascade, the shape of the sucticn 
surface of a blade and the velocity distribution on that surface 
will be assumed and the pressure surface of the adjacent blade will 
be sought. The form in which the initial surface and the velocity 
distribution are assujaed and the construction of an adjacent stream- 
line are shown in figure 1. The orthogonal spacing between a 
givaa streamline and the next streamline constructed from it is 
approximated by a spacing 2 taken perpendicular to the given 
streamline . 

- Basic equations . - The continuity equation for two-dimensional 
flow in a channel may be written as 


w = [yV dL ( 1 ) 


where 

w channel weight flow per unit height 

7 wei^it density 

V velocity 

dL line element perpendicular to direction of flow 

o 

All symbols are redefined in appendix A. 

The condition for irrotatianal flow is expressed as 


S = - CT 


( 2 ) 
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2. The curvature 
such that 


C is assumed to vary as in vortex flow 


Xj+1 _ Xi 
c i+l " c i 


The application of the two methods resulting from these 
assumptions will he discussed in a subsequent section. 

In order to minimize the inaccuracies due to these assumptions, 
they are applied to only a part of the total channel carrying the 
specified fraction n of the total weight flow. Equation (4) can 
then be solved for the unknown velocity parameter Xi+1 and the 

approximate 2 computed from the irrotatiaa-flow equation. With 
this value known, the streamline ^ i+1 can be graphically con- 
structed and its curvature, in turn, used for the construction of the 
next streamline. 

Solution for assumption of constant curvature . - If C is 
assumed equal "to C^, equation (4) becomes 


-Ci nf (w) 


%+l 

Xi 


1 

o k- 1 

(l - x 2 ) ax 


(5) 


The integral can be expressed as 


^ nf (w) = 


«i+l JL 


0 


(1 - X 2 ) 


1_ 

>1 

dX - J 


'Xi A- 


1 

k-1 

(i - x 2 ) ax 


( 6 ) 


■where the unknown is the value X^ + ^ representing the velocity 
on the streamline forming the boundary of the stream, filament carry- 
ing the weight flow specified by nf(w). 
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Far a known value of k, the integral 


•f*r 


g(3r) =J o 


( 1 - X 2 ) 


1 

k-1 

cLX 


can "be computed as a function of the limit This computation 

has "been carried out in the maimer illustrated in appendix B and 
the results tabulated in tables I and H. 

Equation (6) can be written as 


-Cjnf(v) = g(Xi +1 ) - s(X±) (7) 

With the known curvature C ± , the assigned fraction of the wei^it 
flow nf(v), the known velocity parameter 1^, and giX^) found 
from the tabulated values, equation (7) is solved for g(Xi+i) 
and X 1+1 found from table I. 

If X i+1 is known, the approximate spacing l Q (subscript 

denoting constant-curvature assumption) can be found for the 
assumption of constant curvature by integrating equation (2). 

^i+1 X. ,-i 

l0 «e vf = -°i l c ' lo 6o ^ 

or 

1 ^i+1 

l 0 - - 4 («) 

The value Z c is confuted at a point on the known streamline 
and graphically laid out normal to the streamline at that point. 

By computing 2 for a number of points, the adjacent stream- 
line \|/ i+1 can be drawn. (See fig. 1.) With the curvature of 

this new streamline C^ + 2 , which can be readily computed, the 
initial curvature C^, and the spacing 2 , a closer approximation 
to the velocity cn the new streamline 2^ +1 can be found from 
equation (8) using an average curvature. 
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- ( C 1 + C l+l) l 
* *i e 2 


(9) 


Further refineanent of l Is not required unless the average 
curvature Is greatly different from the Initial value C^. If the 

velocity parameter is known, the next streamline 1+2 

he constructed. The process is repeated until channel space has 
heen c omp uted to allow for all the specified weight flow. 

Solution for assumption of vortex variation . - A vortex 
variation of the form — — — 


Ci Ci 


( 10 ) 


Is Inserted Into equation (4). 


r*. 


nf(w) 


1+1 


1 

k-1 


h 1 (i - x*) az 

Cj[ X 


If the constant is removed from the integrand and the 

integral is replaced hy the sum of two Integrals , 

1 1 


- * nf (w) 
*1 


*1+1 FI 

1 (i - x 2 ) ax - 
x 


** FI 

i (l - x 2 ) ax (li) 

A 


At the lower limits the integrals are improper. As only the 
differences between the two integrals are required , a lower limit 
of 0.01 will therefore he used. By definition, 


h(X ) 
r 



i (1 - X 2 ) 
X 


1_ 

k-1 


ax 


.01 
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For a inown value of k, the Integral h(X r ) can he 
confuted as a function of the limit Xp. This computation has 

been performed In the manner Illustrated In appendix C and the 
results tabulated In tables HI and IV. 


Equation (ll) can then be written 

- nf(v) = h(I^ +1 ) - hdj) (12) 


and the value of may be found. 

The channel spacing can then be computed by integrating 
equation (2) with the curvature variation of equation (10) 



or 




By Integration between Xj and X^^ 


_1 (J± . 

Cl ( Xi+i 




(13) 


where the subscript v Indicates vortex-type -variati an. assumption.' 

With this value for l y , the new streamline can be constructed 
and a closer approximation to the velocity parameter can 

be obtained by use of equation (9). 

Comparison of curvature assumptions . - Far the mum Initial 
condition; the effect of the, curvature assumption, on the computed 
value of 2 is a function of the step size so that far a specific 
application the effect is dependent on the number of steps l/n 
used In the computation. Without a knowledge of the true variation 
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In the streamline curvature, the absolute effect of the assumption 
cannot he determined. From several typical design problems, the 
vortex-type variation was found to more closely approximate the 
variations obtained. 

For the same initial conditions of ^ and Xj and for the 

MWfl value of nf(w) , the differences in Z for the constant 
curvature and vortex -type-variation assumptions can be computed. 
These differences are plotted in figure 2 as a function of step- 
size par ame ter C ± Z y with Xj asa parameter. It can be seen 

that the constant-curvature assumption yields a higher value of 
Z than does the vortex-type variation; For most design problems, 
the ma ximum value of C ± l y (or C ± Z c ) can be kept below 0.3 

without an excessive number of steps. For this value, the differ- 
ences in Z are usually less than the accuracy of the graphical 
construction process. 


APPLICATION CEB* METHOD 

Selection of prescribed velocity and boundary shape . - In 
blade -design, work, the suction surface is usually the more critical 
surface. As the proposed method allows any streamline to be 
specified, it appears to be of greater advantage to fix the suction- 
surface conditions. Any arbitrary surface and velocity distribution 
may be assigned, but they might not result in a physically possible 
blade. As a basis for assigning the initial surface and velocity, 
it has been found convenient to analyze first an arbitrary channel 
drawn to give the appearance of a conventional blade. The vel- 
ocities in this channel can be easily found by the direct stream- 
filament method of reference 4. The suction-surface shape or the 
velocity distribution or both can then be adjusted and utilized 
as the initial values in the proposed method. (Such adjustments 
should be caref ully made to prevent excessive thickness or thinness 
of the resulting blade. For subsonic flow, an increase in the 
assigned velocity will result in a narrower channel, thus a thicker 
blade, and vice versa. An increase in curvature on the suction 
surface will result in a wider channel at that point.) 

Boundary-layer theory has shown that flow separation is most 
likely to occur in a region where deceleration takes place. In 
order to minimize the danger of separation, information on boundary- 
layer flow should be employed to limit the rate of flow deceleration 
along the suction surface. In the absence of any reliable infor- 
mation on boundary layers, the velocity over as large a portion of 
the suction surface as possible should be either constant or 
increasing. 
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Determination of weight -flow parameter . - The weight -flow 
P araneter f (v) nay be readily confuted from the known velocity 
diagrams and the blade spacing by 


/ 2 \ k-1 

f(w) = sXoojl - Zee ] COB e (14) 

(See appendix D for development. ) 
w! lere 


s blade pitch 

Z<» relative inlet -velocity parameter at Infinity 

6 angle between relative velocity and Its axial component at 

Infinity 

The parameter f(w) is related to the mass-flow parameter u 
of reference 4 by 



where Iy Is the total chann el width measured a l on g an orthogonal 
line. 


For design problems where the flow per unit blade heigit is 
not constant through the ch ann el, the weight-flow parameter can be 
varied to conform to the expected distribution In weight flow per 
unit blade height. 

Selection of number of steps . - The number of steps to be used 
In the computation can be set from considerations of the desired 
accuracy. As the errors involved In a curvature assumption Increase 
with step size, it is apparent that small steps are desirable: 
as the graphical construction of each streamline introduces errors, 
however, the least number of steps consistent with theoretical 
accuracy should "be used* 

A large number of small steps would, theoretically, yield the 
highest accuracy. Two factors contributed errors to the computation: 
namely, (1) the assumption on the curvature variation and (2) the 
use of a straight line In place of a true orthogonal for construct- 
ing Z. Without knowing the true curvature variation, complete 
©valuation of the effect of these assumptions is impossible. 
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With low initial velocities and low curvatures, relatively 
larger steps can "be used than for higi velocities and curvatures. 

For most turbine blades, five steps (n = l/5) should be adequate; 
whereas thr ee steps can be used for sections not requiring accurately 
specified velocities. 

S election of curvature variation . - The comparison of l as 
computed for the vortex variation and for constant curvature is 
shown in figure 2. For low Initial velocities and curvatures, the 
difference is very small; unless a high degree of accuracy is 
desired and large steps used, the differences are small compared with 
the errors involved in graphically constructing l and computing 
curvatures. From illustrative examples, the vortex variation 
was found to more closely approximate the measured variation than 
did the constant -curvature assumption. 


SUGGESTED PROCEDURE FOR COMPUTATION 
Thft following step-by-step procedure is suggested: 

1. Assign suction-surface shape H'q an ^ 'velocity distribution 
Xq. The curvature of -Hie assigned streamline C Q can be found by 
graphical or computational methods. 

2. Compute the weight-flow function f(v) and set n for 
the desired number of steps. For a departure from two-dimensional 
flow where the value of f(w) is not constant, the values at the 
inlet and the outlet of the channel should be computed and a 
variation between these values assumed. 

3. For a number of points (10 to 15) on the given str eamlin e, 
compute l . 

(a) For vortex variation . - For the initial vel ocity 
parameter Xq, oVfeain the value of 1i(Xq) from table ILL. Ry use 

of equation (12) compute h(X^) and obtain X^ from table 111. 
Compute Z v by equation (13). 

(b) For c ons tant curvature. - For the initial velocity 
parameter Xq, obtain the value of g^) from table I. Ry use 

of equation (7) compute g(X^) and obtain X^ from table X. 

Compute l c by equation (8). 
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4. With the values of 2 known, construct the streamline \J/j 
as shown In figure 1 and compute its curvature C-^. 

5. Compute the velocity an the new stre amlin e by equation (9). 
The curvature and velocity values in this equation are those at the 
points of intersection of the 2 line with the respective stream- 
lines. 


6. Repeat steps 3 and 4 for the assigned number of stream 
filaments and take the last stre amlin e ^ ±/ n as the channel 

boundary. 

7. Compute the velocities in the channel to determine the 
degree to which the specified surface velocities have been satisfied. 
The method of reference 4 can be used to good advantage although it 
may be necessary to carry the method beyond the usual first approxi- 
mation of consideration of the channel in a single step. 


EFFECT OF V-ARXATICNS IN PRIMARY PARAMETERS 


As an aid in adjusting assigned velocity and curvature values 
as well as an indication of the needed accuracy in curvature determi- 
nation, a study was made of the effect of variations in the primary 
parameters Cj[ and an 2 y . 


The effect of a variation in initial curvature on the computation 
of 2 V was found by taking the partial derivatives of equations (12) 

and (13) with respect to C^. As 

bl „ LI 
~ AC a 


and by definition 

LI 

E2 2 
EC* = ACj 

”®i 


it can be shown that 
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Values of this ratio EZ V /EC^ have "been computed far values of 
from 0 to 1.0 and X^ of 0.2, 0.3, and 0.4 for k of 1.4. 

The results are shown in figure 3. 

Similarly, the effect of a variation In Xj an the computation 
of l r was found to he 
AZy 

2i IS 

“ AXj 



Values of this ratio EZy/EX^ have been computed for values of 
CjZ 7 from 0 to 1.0 and Xj of 0.2, 0.3, and 0.4, for k of 1.4. The 
results are shown In figure 4. 

Illustrative Example 

The channel portion of a typical blade was designed by the 
methods presented herein. Assigned values were 
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Xa, = 0.2062 LrSa. = 0.5051 


'cr 


12 = 1.475 
s 


e = 4i. 3 C 


■where 


axial width 


An arbitrary suction-surface velocity and suction-surface shape 
were assumed as X Q and \| /q, as shown in figure 5. For n = l/3, 

the channel was constructed using "both the vortex-type ■variation 
and the constant -curvature method. The design plot for the constant- 
curvature method is shown in figure 5. The third stre amlin e \J/ 3 is 

taken as the pressure surface of the adjacent blade. A comparison 
of the channels found by the two methods is shown in figure 6. In 
the channel section corresponding to the highly curved portion of the 
suction surface, the channel width was greater for the constant- 
curvature computation. This difference is in accordance with the 
"theoretical differences previously noted. Beyond the 40-percent 
chord position, the two channels are in close agreement. The 
differences in pressure-surface profile, which show a wider channel 
for the vortex variation for a small section of the blade surface. 

are possibly due to errors in graphical construction and in curvature 
computation# 


Lewis Flight Peopulsion Laboratory, 

National Advisory Committee for Aeronautics, 
Cleveland, Ohio, March 17, 1949. 
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appmdix a 


STMBOIS 


The following symbols are used In this report: 
curvature (reciprocal of radius of curvature) 
specific heat at constant pressure 
axial width 

prefix to denote fractional variation 


f(v) weight -flow parameter, 


g(X r ) integral function 


7' \/2c^T' 


, /V 


hfXj.) integral function. 


1(1 - X‘ 

x.\ 


ratio of specific heats, Cp /Oy 

length a3 nng velocity, potential line 
approximate width of stream fila m ent 
specified fraction of wei$it flew 
hlade pitch 
absolute temperature 
velocity 

Channel weight flow per unit heigit 
dimensionless velocity parameter, 


n/^V 7 " 
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7 veigbt density 

A prefix to indicate change 

0 angle between relative inlet velocity and its axial 

component at inf ini ty 

^ designation of streamlines 

Subscripts : 

0 assigned conditions 
1,2,3 index of streamlines 

c computed for constant curvature 

\ 

cr conditions at critical sonic velocity 

1 general index 

r limit value 

v computed for vortex-type variation 

w total width 

x axial 

Superscripts : 

’ stagnation state 
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APFMDIX B 

COMPUTATION OF g^) 


The Integral 




ax 


csn "be readily integrated by expanding the integrand In a power 
series flr> a integrating tennwise- 



(2-k)(3-2k) x 6 + 
6(k-l) 3 


g(3^) = 


1 

Ea 

(i - x 2 ) ax = Xp - 


Xy 5 + (2-k)X r b 
3 ( k " 1 ) 10(k-l) 2 


(2-k)(5-2k)X r 7 

42(k-l) b 


Par a tfe.ch number of 1.0 and a k of 1.40, the value of X is 
0.40825 for lower values of k, the value of X decreases. For 

the range of X from 0 to 0.41 and for k from 1.28 to 1.40, the 
first four terms shown In the previous expansion were sufficient to 
compute gCXj.) accurately to five decimal places. 


The values of the integral g(Xy ) are given in table I for 
values of the argument X^. from 0 to X^ in increments of 0.005 
f ar k of 1.40, 1.38, 1.36, 1*34, 1.32, 1.30, and 1.28. 

Values of g(X r ) for critical (sonic) velocities for the values 
of k used are tabulated in table II . For critical flow 


X, 


cr 
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APPENDIX C 


CO^XJEATICN OF hCXj.) 

The integral 

r * r JL 

/ _ k-l 

' MXp) = / I (i - x 2 ) ax 
0.01 x 

can he readily integrated by expanding the integrand in a power 
series and integrating termwise. 


1 

X 


(1 - X 2 ) 


1 

k-l 



1 

k-l 


X + 


2-k y 3 
2(k-l) 2 


(2-k)(5-2k) x 5 + 
6(k-l) 3 


... dX 


When this equation ia integrated, it yields 


h(X r ) 


= log ° 1 - 2irar * 


(2-k) y4 . (2-k)(3-2k) yS J 
8(k-l) 2 36(k-l) 3 


-fr 

0.01 


Values of this integral for X. from 0.010 to X have 

r cr 

■been computed and are given in table HI for k = 1.40, 1.38, 1.36, 
1.34, 1.32, 1.30, and 1.28. 

Values of h(X^) for critical velocities are given in table IV # 
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APPENDIX D 


DKEEHa-HHATICN OF WEIGH!! -FLOW PARAMETER f (w) 

Thft weight -flow parameter may "be readily determined from the- 
velocity diagrams by application of the continuity equation to can- 
•ditians at infinity relative to the blade. The channel weight flow 
per unit blade height may be written as 


w 


ys V. 


X,co 


The term f(w) is defined as 


f(w) = 


w 

7 * ^cJF 


so that 


Inasmuch as 




f ( w ) = JL -1^12 


and 


JL = 

7 ' 


1 


k-1 f y\ 2 

" k+1 XycrJ 


1 

FI 


V v = Yoc cos e 

AyCD 


and 


X = 

si k+1 Y cr 
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then 


f (w) = sXjl - x c 


1 

lk-1 


cos 6 


For rotor "blades f (w) must be computed relative to the 
rotating blades. 


REFERENCES 

1. Mutterperl, William: A Solution of the Direct and Inverse 

Potential Problems for Arbitrary Cascades of Airfoils. 

NACA AKR L4K22b, 1944. 

2. Goldstein, Arthur W. , and Jerison, Meyer: Isolated and Grid 

Airfoils with Prescribed Velocity Distribution. NACA TN 1308, 


3. Stodola, A. : Steam and Gas Turbines. Vol. H. McGraw-Hill 

Book Co., Inc. , 1927, p. 992. (Reprinted, Peter Smith (New 
York), 1945.) 

4. Huppert, M. C., and MacGregor, Charles: Comparison between 

Predicted and Observed Performance of Gas-Turbine Stator Blade 
Designed for Free-Vortex Flow. NACA TN 1810, 1949. 

5. von Mises, Richard, and Friedrichs, Kurt 0. : Fluid Dynamics. 

Advanced In s true t i on and Research in Mechanics, Brown Uhiv* , 
Summer 1941, p. 30. 

6. Weinig, F. : Flow Patterns of a Compressible Gas through Turbine 

Cascades Whoa Shock Waves are Avoided. Tech. Rep. 

No. F-TR-2141-ND, Air Materiel Command, Army Air Forces, 

May 1947. 



22 


NACA TN 1931 



.00500 

.01000 

.01500 

.01999 

.02499 

.02998 

.03496 

.03995 

.04492 

.04990 

.05486 

.05982 

.06477 

.06972 

.07465 

.07958 

.08449 

.08940 

.09429 

.09917 

.10404 

.10890 

.11374 

.31857 

.12338 

.12818 

.13297 

.13773 

.14248 

.14722 


1.38 

1.36 

1.34 

1.32 

1.30 

1.28 

0 

0 

0 

0 

0 

0 

.00500 

.00500 

.00500 

.00500 

.00500 

.00500 

.01000 

.01000 

.01000 

.01000 

.01000 

.01000 

.01500 

.01500 

.01500 

.01500 

.01500 

.01500 

.01999 

.01999 

. .01999 

.01999 

.01999 

.01999 

.02499 

.02499 

.02498 

.02498 

.02498 

.02498 

.02998 

.02998 

.02997 

.02997 

.02997 

.02997 

.03496 

.03496 

.03496 

.03496 

.03495 

.03495 

.03994 

.03994 

.03994 

.03993 

.03993 

.03993 

.04492 

.04492 

.04491 

.04490 

.04490 

.04489 

.04989 

.04988 

.04988 

.04987 

.04986 

.04985 

.05485 

.05485 

.05484 

.05483 

.05482 

.05480 

.05981 

.05980 

.05979 

.05978 

.05976 

.05974 

.06476 

.06475 

.06473 

.06472 

.06470 

.06467 

.06970 

.06968 

.06966 

.06964 

.06962 

.06960 

. 07463 > 

.07461 

.07459 

.07456 

.07453 

.07450 

.07955 

.07953 

.07950 

.07947 

.07943 

.07939 

.08446 

.08443 

.08440 

.08436 

.08432 

.08427 

.08936 

.08933 

.08929 

.08924 

.08920 

.08914 

.09425 

.09421 

.09416 

.09411 

.09405 

.09399 

.09913 

.09908 

.09902 

.09896 

.09890 

.09882 

.10399 

.10393 

.10387 

.10380 

.10372 

.10363 

.10884 

.10878 

.10870 

.10862 

. 10853 

.10843 

.11367 

.11360 

.11352 

.11343 

.11333 

.11321 

.H850 

.11841 

.11832 

.11822 

.11810 

.11797 

.12330 

.12321 

.12310 

.12299 

.12285 

.12270 

.12809 

.12798 

.12787 

.12774 

.12759 

.12742 

.13286 

.13274 

.13261 

.13247 

.13230 

.13211 

.13762 

.13749 

.13734 

.13718 

.13699 

.13678 

.14235 

.14221 

.14205 

.14187 

.14166 

.14143 

.14707 

.14691 

.14673 

.14653 

.14631 

1 .14605 
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Eatio of specific heats of gas, k 


1.38 1.36 1.34 1.32 1.30 




.16952 

.17405 

.17855 

.18302 

.18746 



1.28 


4653 0.14631 0.14605 
5118 .15093 .15065 
5580 .15553 .15522 

.16040 .16010 .15976 
.16498 .16465 .16428 


.16917 

.17366 

.17813 

.18257 

.18698 

.19136 
.19571 
.20002 
S' .20431 
.20856 

.21278 

.21697 

.22113 

.22525 

.22933 

.23338 

.23740 

.24138 

.24532 

.24922 

.25309 

.25691 

.26070 

.26445 

.26816 

.27184 


.16877 

.17323 

.17766 

.18206 

.18643 

.19077 

.19507 

.19934 

.20358 

.20779 

.21196 

.21610 

.22020 

.22426 

.22828 

.23226 

.23623 

.24014 

.24401 

.24785 

.25164 

.25540 

.25912 

.26279 

.26642 

.27001 
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X r 1 

r 2 "ETI 

TABLE I.- VALUES OF FUNCTION g(X r ) - I (1-X ) dX - Concluded, 

0 


Velocity 


Ratio of specific heats of gas, 

k 


pCUTcuH© U ©r 

ly . 

1.40 

1.38 

1.36 

1.34 

1.32 

1.30 

1.28 • 

0.300 

0.27840 

0.27734 

0.27618 

0.27489 

0.27345 

0.27184 

0.27001 

.305 

.28234 

.28123 

.28001 

.27866 

.27715 

.27547 

.27356 

.310 

.28624 

.28508 

.28380 

.28239 

.28082 

.27906 

.27707 

.315 

.29010 

.28889 

.28756 

.28609 

.28445 

.28261 

.28053 

.320 

.29394 

.29268 

.29128 

.28975 

.28803 

.28612 

.28395 

.325 

.29774 

.29642 

.29497 

.29337 

.29158 

.28958 

.28733 

.330 

.30150 

.30013 

.29862 

.29695 

.29509 

.29301 

.29066 

.335 

.30523 

.30380 

.30223 

.30049 

.29856 

.29639 

.29396 

.346 

.30893 

.30744 

.30580 

.30399 

.30198 

.29973 

.29720 

.345 

.31259 

.31104 

.30934 

.30746 

.30537 

.30303 

.30040 

.350 

.31621 

.31460 

.31284 

.31088 

.30871 

.30629 

.30356 

.355 

.31980 

.31813 

.31629 

.31427 

.31202 

.30950 

. .30667 

.360 

.32335 

.32162 

.31971 

.31761 

.31528 

.31267 


.365 

.32687 

.32507 

•3231C 

.32092 

.31850 

.31580 


.370 

.33035 

.32848 

.32644 

.32418 

.32167 



. 375 

.33379 

•33186 

.32974 

.32740 

.32481 



.380 

.33719 

.33520 

.33300 

.33058 




.385 

.34056 

.33849 

.33622 

.33372 




.390 

.34389 

.34175 

.33940 





.395 

.34718 

.34497 

.34255 





.400 

.35043 

.34815 





1 

.405 

.35365 







.410 

.35682 








TABLE II.- VALUES OF g(X r ) FOR CRITICAL VELOCITIES 


k 

*cr 

s(x cr ) 

1.40 

1.38 

1.36 

1.34 

1.32 

1.30 

1.28 

0.40825 

.39957 

.39056 

.38118 

.37139 

.36115 

.35044 

0.35571 

.34788 

.33976 

.33132 

.32254 

.31339 

.30384 
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Xr l 

TABLE III - VALUES OF FUHCTION h(Xy) • J" iU-X 2 )^ dX 


O.Ol 


Velocity 

parameter 

Ratio of specific heats of gas, k 

- h . 

1.40 

1.38 

1.36 

1.34 

1.32 

1.30 

1.28 

0.010 

0 

0 

0 

0 

0 

0 

0 

.015 

.40531 

.40530 

.40529 

.40528 

.40527 

.40526 

.40524 

.020 

.69277 

.69275 

.69273 

.69271 

.69268 

.69265 

.69261 

.025 

.91564 

.91600 

.91556 

.91552 

.91547 

.91542 

.91535 

.030 

1.09761 

1.09756 

1.09750 

1.09744 

1.09736 

1.09728 

1.09718 

.035 

1.25136 

1.25128 

1.25120 

1.25111 

1.25101 

1.25089 

1.25076 

.olio 

1*38442 

1.38432 

1.38421 

1.38409 

1.38395 

1.38380 

1.36362 

.045 

1.50167 

1.50155 

1.50141 

1.50125 

1.50107 

1.50087 

1.50064 

.050 

1.60644 

1.60628 

1.60611 

1.60591 

1.60569 

1.60544 

1.60516 

.055 

1.70110 

1.70090 

1.70069 

1.70045 

1.70018 

1.69988 

1.69954 

. 060 

1.78739 

1.78716 

1.78691 

1.78662 

1.78630 

1.78594 

1.78552 

.065 

1.86665 

1.86638 

1.86608 

1.86575 

1.86537 

1.86494 

1.86446 

.070 

1.93992 

1.93961 

1.93926 

1.93887 

1.93843 

1.93793 

1.93737 

.075 

2.00801 

2.00765 

2.00725 

2.00680 

2.00630 

2.00572 

2.60567 

.080 

2.07159 

2.07117 

2.07072 

2.07021 

2.06963 

2.06898 

2.06824 

.085 

2.13118 

2.13072 

2.13020 

2.12962 

2.12898 

2.12824 

2.12740 

.090 

2.18726 

2.18673 

2.I8615 

2.18551 

2.18478 

2.18396 

2.18301 

.095 

2.24017 

2.23959 

2.23895 

2.23822 

2.23741 

2.23650 

2.23545 

.100 

2.29026 

2.28961 

2.28890 

2.28810 

2.28720 

2.28618 

2.28502 

.105 

2.33778 

2.33706 

2.33628 

2.33540 

2.33441 

2.33328 

2.33200 

.110 

2.38296 

2.38218 

2.38132 

2.38035 

2.37927 

2.37804 

2.37663 

.115 

2.42602 

2.42517 

2.42423 

2.42317 

2.42198 

2.42064 

2.41911 

.120 

2.46713 

2.46620 

2.46517 

2.46403 

2.46274 

2.46127 

2.45961 

.125 

2.50644 

2.50543 

2.50432 

2.50307 

2.50167 

2.50009 

2.49828 

.130 

2.54408 

2.54300 

2.54179 

2.54045 

2.53894 

2.53723 

2.53528 

.135 

2.58019 

2.57902 

2.57772 

2.57627 

2.57464 

2.57280 

2.57070 

,l4o 

2.61486 

2.61360 

2.61221 

2.61065 

2.60891 

2.60693 

2.60467 

.145 

2.64820 

2.64685 

2.64536 

2.64369 

2.64182 

2.63970 

2.63729 

.150 

2.68029 

2.67885 

2.67725 

2.67547 

2.67347 

2.67121 

2.66863 


■ — v — 
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TABLE ITT - VALUES OF FUHCTICB bClj.) 


^X_ ' 1 

/ 

O.Ol 


dX - Continued. 


Velocity 

na ra mot.PY* 


Batlo of speclfio heats of gas 

4 


x r 

1.40 

1.38 

1.36 

1.34 

1.32 

1.30 

1.28 

0.150 

2.68029 

2.67885 

2.67725 

2.67547 

2.67347 

2.67121 

2.66863 

.155 

2.71120 

2.70967 

2.70797 

2.70607 

2.70393 

2.70152 

2.69877 

.160 

2.74102 

2.73939 

2.73758 

2.73554 

2.73329 

2.73072 

2.72780 

.165 

2.76980 

2.76807 

2.76614 

2.76400 

2.76159 

2.75887 

2.75577 

.17© 

2.79760 

2.79576 

2.79373 

2.79145 

2.78890 

2.78602 

2.78273 

.175 

2.82448 

2.82254 

2.82038 

2.81798 

2.81528 

2.81223 

2.80876 

.100 

2.85049 

2.84843 

2.84616 

2.84362 

2.84077 

2.83755 

2.83388 

.185 

2.87566 

2.07350 

2.87109 

2.86842 

2.86541 

2.86202 

2.85816 

.190 

2.90005 

2.89777 

2.89524 

2.89242 

2.88926 

2.88569 

2.88163 

.195 

2.92368 

2.92128 

2.91863 

2.91567 

2.91235 

2.90860 

2.90433 

.200 

2.94660 

2.94409 

2.94130 

2.93819 

2.93470 

2.93077 

2.92629 

.205 

2.96884 

2.96620 

2.96328 

2.96002 

2.95637 

2.95225 

2.94756 

.210 

2.99043 

2.98766 

2.98460 

2.98119 

2.97737 

2.97306 

2.96815 

.215 

3.01139 

3.00850 

3.00530 

3.00173 

2.99774 

2.99323 

2.98810 

.220 

3.03176 

3.02874 

3.02539 

3.02167 

3.01750 

3.01279 

3.00744 

.225 

3.05155 

3.04840 

3.04491 

3.04102 

3.03667 

3.03176 

3.02618 

.230 

3.07080 

3.06751 

3.06387 

3.05982 

3.05529 

3.05017 

3.04436 

.235 

3.08952 

3.08608 

3.08230 

3.07809 

3.07336 

3.06804 

3.06200 

.240 

3.10772 

3.10416 

3.10022 

3.09583 

3.09092 

3.08539 

3.07911 

.245 

3.12544 

3.12174 

3.11764 

3.11309 

3.10799 

3.10224 

3.09571 

.250 

3.14270 

3.13885 

3.13459 

3.12986 

3.12457 

3.11860 

3.11183 

.255 

3.15949 

3.15550 

3.15109 

3.14618 

3.14069 

3.13450 

3.12748 

.260 

3.17585 

3.17171 

3.16714 

3.I6205 

3-15636 

3.14995 

3.14268 

.265 

3.19178 

3.18750 

3.18276 

3.17749 

3.17160 

3.16496 

3.15744 

.270 

3.20731 

3.20287 

3.19796 

3.19251 

3.18642 

3.17956 

3.17178 

.275 

3.22244 

3.21785 

3.21277 

3.20714 

3.20084 

3.19375 

3.18571 

.280 

3.23719 

3.23244 

3.22720 

3.22137 

3.21486 

3.20754 

3.19925 

.285 

3.25156 

3.24666 

3.24125 

3.23523 

3.22851 

3.22096 

3.21240 

.290 

3.26558 

3.26052 

3.25493 

3.24872 

3.24179 

3.23400 

3.22518 

.295 

3.27925 

3.27403 

3.26827 

3.26187 

3.25472 

3.24669 

3.23760 

.300 

3.29258 

3.28720 

3.28126 

3.27466 

3.26730 

3.25903 

3.24967 
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*r 1 

TABLE III - VALUES OF FUBCTIC* h(l r ) = j Id-X 2 ) 1 " 1 dl . Concluded. 

O.Ol 


Velocity 

parameter 


Ratio of specific beats of gas 

3 , k 


*r 

l.to 

1.38 

1.36 

1.3* 

1.32 

1.30 

1.28 

0.300 

3.29258 

3.28720 

3.28126 

3.27*66 

3.26730 

3.25903 

3.2*967 

.305 

3.30558 

3.3000* 

3.29392 

3.28713 

3.27955 

3.27103 

3 . 261*0 

.310 

3.31827 

3.31256 

3.30626 

3.29927 

3.291*7 

3.28271 

3.27280 

•315 

3.33065 

3.32477 

3.31829 

3.31110 

3.30308 

3.29*07 

3.28389 

.380 

3.3*272 

3.33668 

3.33001 

3.32262 

3.31*37 

3.30512 

3.29*66 

.385 

3.35*50 

3.3*829 

3.3*14* 

3.3338* 

3.3*537 

3.31587 

3.3051* 

.330 

3.36600 

3.35962 

3.35258 

3.34*78 

3.33608 

3.32633 

3.31532 

.335 

3.37722 

3.37067 

3.363*4 

3.355*3 

3.3*651 

3.33651 

3.32521 

.3*0 

3.38817 

3.38144 

3.37*03 

3.365a 

3.35666 

3.3*641 

3.33*83 

.3*5 

3.39886 

3.39196 

3.38*35 

3.37593 

3.36655 

3.35604 

3.3*418 

.350 

3 . *0929 

3.40221 

3.39**2 

3.38578 

3.37617 

3.365*1 

3.35327 

.355 

3. *19*7 

3.41222 

3.40*23 

3.39538 

3.38554 

3.37*52 

3.36210 

.360 

3. *29*0 

3.42198 

3.41379 

3.40*7* 

3.39467 

3.38339 


.365 

3. *3910 

3.43150 

3.42312 

3.41385 

3.40355 

3.39201 


.370 

3. **857 

3.4*078 

3.43221 

3.42273 

3.41219 



.375 

3.45781 

3.4*98* 

3.4*108 

3.43138 

3.42061 



.380 

3.46683 

3.45868 

3-44972 

3.43981 




.385 

3.47563 

3.46730 

3.4581* 

3.44802 




.390 

3.48422 

3.47571 

3 .46636 





.395 

3.*926l 

3.48391 

3.47*36 





.*00 

3.50079 

3.49192 






.*05 

3.50878 







.*10 

3.51657 








TABLE IV - VALUES OF h(l r ) FOE CRITICAL -VELOCITIES 


k 

X„„ 

cr 

*U C r) 

1.40 

1.38 

1.36 

1.3* 

1.32 

1.30 

1.28 

0.40825 

.39957 

.39056 

.38118 

.37139 

.36115 

.350*4 

3.51384 

3.49123 

3.46725 

3.4*175 

3.41453 

3.38537 

3.35*04 
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Projected chord, percent 


Figure I. - Typical assigned suction-surface contour 
typical assigned velocity distribution, and 
construction of first stream filament. 


Channel width for constan t curvature 
Channel width for vortex-type variation 



Step-size parameter, C ; i v 

Figure 2, - Effect of step-size parameter on ratio 
of channel width for constant curvature to channe 
width for vortex-type variation. Ratio of specil 
heats k f 1 . 4 , 


30 
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Figure 3. - Effect of step-size parameter on ratio 

of percentage variation in channel width for 
vortex-type variation to percentage variation in 
initial curvature from equation (15). Ratio of 
specific heats k, 1 . 4 . 



Percentage variation in channel width for vortex-type va r i a t i o n El 
Percentage variation in initial velocity parameter * EX- 
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Step-size parameter, C } 1 v 



Figure 4. - Effect of step-size parameter on ratio 
of percentage variation in channel width for 
vortex-type variation to percentage variation in 
initial velocity parameter from equation (16). 
Ratio of specific heats k, 1 . 4 . 



Projected chord, percent 

Figure 5. - Assigned suet i on -surface contour, assigned velocity distribution 
and channel construction for illustrative example. 









